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The observation of twenty cosmic-ray air-showers at and above 10 eV poses fascinating problems for 
particle astrophysics: how the primary particles are accelerated to these energies, how the primaries 
get here through the 2.7K microwave background filling the Universe, and how the highest-energy 
events exhibit clustering on few-degree angular scales on the sky when charged particles are expected 
be bent by cosmic magnetic fields. An overview of the puzzles is presented, followed by a brief 
discussion of many of the models proposed to solve these puzzles. Emphasis is placed on (i) the 
signatures by which cosmic ray experiments in the near future will discriminate among the many 
proposed models, and (ii) the role neutrino primaries may play in resolving the observational issues. 
It is an exciting prospect that highest-energy cosmic rays may have already presented us with new 
physics not accessible in terrestrial accelerator searches. 
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I. INTRODUCTION 

An unsolved astrophysical mystery, now forty years 
old, is the origin and nature of the extreme energy cos- 
mic ray primaries (EECRs) responsible for the observed 
events at highest energies, 10^" eV [1]. About twenty 
events at ~ 10^" eV have been observed by five different 
experiments [2] . The origin of these events is a mystery, 
for there are no visible source candiates within 50 Mpc 
except possibly M87, a radio-loud AGN at ~ 20 Mpc, 
and Cen-A (NGC5128), a radio galaxy at 3.4 Mpc, and 
neither of these is in the direction of any of the observed 
events [3] . Since the observed events display a large-scale 
isotropy, many sources rather than one source seem to be 
required. The nature of the primary particle is also mys- 
terious, because interactions with the 2.73K cosmic mi- 
crowave background (CMB) renders the Universe opaque 
to nucleons at Eqzk ^ 5 x 10^^ eV, and double pair pro- 
duction on the cosmic radio background (CRB) renders 
the Universe opaque to photons at even lower energies. 
The theoretical prediction of the end of transparency for 
nucleons at Eqzk ~ 0.5 x 10^° eV is the famous "GZK 
cutoff" [4]. Figure 1 shows a recent compilation of the 
AGASA data set, clearly extending beyond Eqzk- 

The main theory challenges in attempting to under- 
stand the super-GZK data are: (i) What cosmic source 
could have accelerated the primary particles to such ex- 
tremely high energies? and (ii) If the sources are distant 

(~ 100 Mpc), then how could their primaries have prop- 
agated through the cosmic background radiation without 
substantial energy loss? The acceleration mechanism ci- 
ther requires a Zevatron accelerator {IZeV = 10^^ eV) 
[5], distant because such a source could not be missed 
if it were nearby; or speculative decaying super-massive 
particles (SMPs) or topological defects (TDs) with mass- 
scale ~ 10^^ eV, clustered nearby; or possibly magnetic 
monopoles accelerated by the cosmic magnetic fields. 
From distant Zevatrons, only neutrinos among the known 
particles can propagate unimpeded to earth. Exploit- 
ing this fact are the "Z-burst model," and the strongly- 



interacting neutrino model. Adding to the drama and 
mystery at present is the observed large-scale isotropy 
and small-scale anisotropy. Surely these characteris- 
tics hint at a solution to the mystery of origin; protons, 
nuclei, and magnetic charges bend in cosmic magnetic 
fields, whereas photons and neutrinos do not. With more 
data above Eqzk, several distinct telltale signatures in- 
cluding the isotropy/anisotropy will allow one to discrim- 
inate among the many models proposed for the origin and 
nature of the EECRs. 
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FIG. 1. Extreme-energy cosmic ray spectrum as observed 
by AGASA. Error bars correspond to 68 % C.L. and the num- 
bers count the events per energy bin. The dashed line reveal- 
ing the GZK cutoff is the spectrum expected from uniformly 
distributed astrophysical sources (from the AKENO website 
[12]). 
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II. EVENT CLUSTERING 

One revealing signature is already evident in the ex- 
isting data sample. This is the pairing of events on the 
celestial sky. The AGASA experiment has already pre- 
sented data strongly suggesting that directional pairing 
is occuring at higher than chance coincidence [6] . Of the 
47 pubUshed AGASA events above 4 x 10^^ eV, 9 are con- 
tained in three doublets and one triplet with separation 
angle less than the angular resolution of 2.5°. The chance 
probability of this clustering occuring in an isotropic dis- 
tribution is less than 1%. The chance probability for 
the triplet alone is only 5%. Of the seven events above 
10^" eV, three are counted among the doublet events. 
Most recently, AGASA has reported two more events 
above Eqzk [7] . Each aligns in direction with a previous 
event, reducing further the probability for random clus- 
tering to 0.07% at 3(7. Comparisons of event directions in 
a combined data sample of four experiments further sup- 
ports non-chance coincidences, especially in the direction 
of the SuperGlactic Plane [8] . 

Whether the pairing is random or dynamical, we shall 
also know in the near future. If the pairing turns out to 
be dynamical, 1 would argue that neutrinos arc a favorite 
candidate for the primary particles. This is because pho- 
tons have such a short 10 Mpc) abosrption length, 
and protons are bent by cosmic magnetic fields during 
their extragalactic journey. 

If neutrinos are the primaries, they should point back 
to their sources, thereby enabling point-source astron- 
omy for the most energetic sources of flux at and above 
10^" eV. It was reported [9] that the first five events at 
10^" eV did in fact point toward extragalactic compact 
radio-loud quasars, just the class of objects which could 
accelerate EECRs to ZcV energies via shock mechanisms. 
With the inclusion of subsequent data, this association is 
controversial [10] . The jury awaits further evidence. 



A. Coincidence or physics? 

A random distribution is obtained by tossing n events 
randomly into ~ {n/Tr0)^ = 1045 (f7/1.0 sr)(6'/1.0°)-2 
angular bins, where fl is the solid angle on the celestial 
sphere covered by the experiment and is the bin half- 
angle. Each resulting event distribution is specified by 
the partition of the n total events into a number mo of 
empty bins, a number mi of single hits, etc., among the 
N angular bins. The probability to obtain a given event 
topology is [11] 



P{{mi},n,N) = 
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The variables in the probability are not all independent, 
as ^j^i j X nij = n and J2j=o ~ ^ useful to 

use these constraints to rewrite this probability as 
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When TV > n > 1, a hmit valid for the AGASA, 
HiRes, Auger and Telescope Array experiments, one finds 
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where r = {N — mo)/n « 1. The non-Poisson nature of 
eq. (4) is reflected in the factorials and powers of r in the 

exponents. 

Typical values of effective area A (km^ sr), celestial 
solid angle fl, and angular resolution ^min for the existing 

and proposed EECR experiments [12] are shown in the 
following Table, where the incident flux F{> Eqzk) = 



10- 
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sr has been used to estimate the number 



of events (n/yr) above Eqzk = 5 x 10^^ eV. 





AGASA 


HiRes 


Auger/TA 


EUSO/OWL 


km^ sr 


150 


800 


6 X lO'^ 


3 X lO'^ 


n/yr 


5 


30 


200 


10* 


n (sr) 


4.8 


7.3 


4.8 


4tt 


^min 


3.0° 


0.5° 


1.0° 


1.0° 



In Figure 2 are shown inclusive probabilities for the 
Auger experiment with 100 events, as determined by for- 
mulae eq. (4). Note that "inclusive probability" means 
the stated number of j-plets plus any other clusters. The 
8-doublet probability is extremely sensitive to the angu- 
lar binning; observation of a flatter dependence on an- 
gular bin-size could signal a non-chance origin for the 
clustering. The observation of two triplets with angular 
binning of less than 2° would constitute 3-sigma evidence 
for cosmic dynamics. The same conclusion would hold if 
a quadruplet within 2.5° is observed among the first 100 
Auger events. 
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Figure 2: Inclusive probabilities for various clusters in a 100 
event sample at Auger. The solid line is the exact result, the 
dashed line is the Poisson approximation 



Without some modification, our analytic formula may 

not be directly relevant to HiRes data. Clear, moonless 
nights are required for detection of atmospheric fluores- 
cence, and summer nights 18 hrs right ascension) are 
efi^ectively 40% shorter than winter nights (~ 6 hrs right 
ascension) for monocular HiRes [13]. Accordingly, the 
HiRes efficiency versus Galactic longitude varies signifi- 
cantly, roughly as N{RA) = N{1 -f e sin(27ri?^/24hrs)), 
with e ~ 0.25. One must ask whether this sinusoidally- 
varying efficiency invalidates the analytic approach with 
its assumption of a constant efficiency. In the Monte 
Carlo approach there is an easy method to generate 
"background" data sets with the experimental efficiency 
properly incuded - one randomly permutes the RA coor- 
dinates of the real data to remove any dynamical correla- 
tion among the events. Fortunately, there is also a simple 
method to estimate the efficiency correction to an ana- 
lytic Poisson distribution, if one assumes that although 
N and n vary with right ascension, ^ = n/N is constant. 
Then it is easy to show that the Poisson distribution 



P{j,m) = 
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For — ^ <^ the exponential is near miity and the 
correction has a closed form, ~ (1— e-^)™/^P^( ^^^^^ ), 
with Pm a Legendre polynomial. As a series, the correc- 
tion is 
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In particular, for <C N, we have Ji = 1 independent 

of e, Xa = 1 + eV2, = 1.03 for HiRes, I3 = 1 + 3eV2, = 
1.09 for HiRes, and I4 = 1 + 3e^ + Se^/8,= 1.189 for 
HiRes. The lesson is that HiRes efficiency corrections 
appear negligible for small numbers of clusters, and we 
may proceed with our analytic analysis for the HiRes 
experiment, for which about 20 events at 10^° eV are 
expected when the first full year's data is analyzed. We 
display in Figure (3) the inclusive probabilities for one or 
more, two or more, and three or more doublets; and one 
or more triplets, over a range of angular binning. 
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Figure 3: Inclusive probabilities for various clusters, given 20 
events at HiRes. The solid line is the exact result, the dashed 
line is the Poisson approximation. 

Note that for all except the 3 doublet configuration, the 
Poisson approximation using the mean values in Eq.(3) 
provides an estimate good to within 50% of the non- 
approximate form; for the (much suppressed) 3 doublet 
configuration, it overestimates the probability by about 
a factor of 3 in much of the angular region. For angular 
binning tighter than 2°. an observation of two doublets 
among the first 20 events has a chance probability of less 
than 0.5%. Thus the observation of this topology could 
be construed as evidence (at the 3a level) for clustering 
beyond statistical. The observation of a triplet within 
< 3° has a random probability of less than 10^'^, and 
hence observation of such a triplet would most likely sig- 
nify clustered or repeating sources, or magnetic focusing 
effects. With the accumulation of 40 events (not shown 
in the figure) , the appearance of two doublets has a prob- 
ability of less than 0.5% for a correlation angle of 1° or 
less. This illustrates how the good angular resolution of 
HiRes may be used to detect non-statistical clustering 
with only a few observed clusters. 

Projected event rates for the EUSO/OWL/AW ex- 
periments present a pleasant problem for our analytic 
formula. In the case where n > iV ^ 1, relevant for 
the EUSO/OWL/AW experiments after a year or more 
of running, higher j-plets are common and the distri- 
bution of clusters can be rather broad in j. From eq. 
(3) we have {rnj/nij-i) = {n/jN) ~ {■nn0'^/j^). Al- 
ready at j = 1 (2), Stirling's approximation to j\ is 
good to 8% (4%), and so for j > 1 we may approxi- 
mate Wj « {N^ /2iTen)^ {en/jNy^^. Extremizing this 
expression with respect to j, one learns that The most 
populated j-plet occurs near j ^ n/N. Combining this 
result with the broad distribution expected for large n/N, 
one expects clusters with j up to several x ^ to be com- 
mon in the EUSO/OWL/AW experiments. Probabilities 
and meaning for frequent, large clusters are somewhat 
difficult to assess numerically. 
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The random distributions displayed in figs. (2) and 
(3) are expected from some models, such as randomly- 
situated decaying super-massive particles (SMPs), or 
charged-particle or magnetic monopole primaries with 
directions randomized by incoherent cosmic magnetic 
fields. A complementary approach to the random prob- 
abilities shown here is to consider specific source mod- 
els generating non-random angular distributions. Steps 
along this line of inquiry have recently been taken [14]. 
Future progress in the field will involve comparisons of 
the random and non-random model predictions with the 
data. 



III. MODELS FOR SUPER-GZK EVENTS 

The conjectured origins of the supcr-GZK particles fall 
into four basic categories. These arc (i) nearby acceler- 
ators, (ii) exotic primaries, (iii) exotic physical law, and 
(iv) neutrino primaries. Other reviews [1] have empha- 
sized the first three categories (especially (i) and (ii)), so 
here we will be brief with those, and put more emphasis 
on the neutrino-primary option. 

A. Nearby accelerators 

Several types of sources have been proposed to exist 
within our Galactic halo. These include highly-ionized 
relativistic dust grains, Galactic supershocks, young neu- 
tron stars, magnetars (highly magnetized pulsars), de- 
caying SMPs with GUT masses of order 10^^ GeV or with 
inflation-motivated masses of order 10^'^^^ GeV, topo- 
logical defects [15] such as strings, Q-balls and vortons, 
and annihilating monopole-antimonopole bound states 
(monopolonium) [16]. For rare sources emitting charged 
particles, such as magnetars, it is necessary to postulate 
that the primaries are iron nuclei to ease the accelera- 
tion requirement and to isotropize the flux in our galactic 
magnetic fleld. It is intriguing that the observed magni- 
tude of CMB fluctuations fixes the reheat temperature 
following inflation to lO^''*^ GeV, which allows gravita- 
tional and thermal production of TDS and/or SMPs of 
just this mass (if they exist) [17]; this mass scale is just 
right for producing 10^° secondary particles via decay. 

For some TD models, dimensional arguments and scal- 
ing laws seem to give an emission rate to short-lived 
SMPs consistent with the observed EECR rate without 
tuning exotic parameters. However, for most decaying 
sources such as SMPs, TDs, and monopolonium, it is 
necessary to tune the lifetime to be longer but not too 
much longer than the age of the Universe in order to 
maintain an appreciable secondary particle emission rate 
today. Discrete gauged symmetries [18] or hidden sectors 
[19] are introduced to stabilize the heavy particle. Then 
rather esoteric physics is needed to break the new sym- 
metry super-softly to maintain the long lifetime. High- 



dimension operators, wormholes [20], and instantons [21] 
bibiteminstanton have been invoked for this purpose. 

Models with sources mainly in the Galactic plane are 
disfavored by the lack of any observed planar asisotropy 
in the data. Models with sources clustered in the Galactic 
halo predict a dipole enhancement in the direction of the 
Galactic center [22], which will be tested by the Auger 
Observatory in a few years. 

Possible sources outside our halo, but still relatively 
nearby, include the radio- loud quasar M87 at 18 Mpc [23] , 
the similar Cen-A at 3.4 Mpc [24], rare nearby GRBs, 
and now-dormant rare AGNs (also called massive dark 
objects - MDOs). These sources are few at best, and 
strong magnetic flelds must be postulated to isotropize 
and/or confine their emissions. In the case of GRBs, 
the identification of the red-shift of their host galaxies at 
typically ^ 1 renders a local occurence of a GRB highly 
improbable [25]. 

Finally, primordial black holes (PBHs) have also been 
suggested, but any sensible initial mass spectrum is un- 
able to provide a sufficient number of PBHs in the final 
stage of decay today. 

B. Exotic primaries 

The GZK cutoff can be raised by simply postulating 
a primary hadron slightly heavier than the proton. The 
reason is kinematical - the cutoff energy varies as the 
square of the mass of the first excited resonant state. 
For this reason, Farrar has proposed that light super- 
symmetric baryons, made from a light gluino plus the 
usual quarks and gluons, may be the primary EECRs 
[26]. Such a scenario renders the source-energetics issue 
even more challenging. In any event, terrestrial exper- 
iments seem to have recently closed the window on a 
possible light gluino. 

Another interesting possibility for the primary EECR 
is the magnetic monopole [27,28]. Any breaking of 
a semisimple gauge symmetry which occurs after in- 
flation and which leaves unbroken a (7(1) symmetry 
group may produce an observable abundance of magnetic 
monopoles. The monopole mass is expected to be ^ a^^ 
times the temperature Tc of the symmetry breaking. At 
the time of the phase transition, roughly one monopole 
or antimonopole is produced per correlated volume. The 
resulting monopole number density today is 

UM ~ 10-^^ (Te/10"GeV)3(ZH/ec)' cm'^, (8) 

where is the phase transition correlation length, 
bounded from above by the horizon size Ih at the time 
of the transition. In a second order or weakly first order 
phase transition, the correlation length is comparable to 
the horizon size. In a strongly first order transition, the 
correlation length is considerably smaller than the hori- 
zon size. 
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The kinetic energy imparted to a magnetic monopole 
on traversing a magnetic field B is Ek ~ gB^, where 
g = e/2a is the magnetic charge according to the Dirac 
quantization condition, and ^ specifies the field's coher- 
ence length. Given the magnitude and coherence length 
data for the cosmic magnetic fields, monopole kinetic en- 
ergies in the range fO^" to 10^"^ cV arc expected; the 
acceleration problem is naturally solved. Monopoles 
with M ~ 10^^ GeV should be relativistic, and carry 
the appropriate energy to qualify as candidates for the 
EECR primaries [28,29]. Within field theory there exist 
many possibilities for an intermediate unification scale 
and intermediate-mass monopoles. 

The propagation problem is also naturally solved. The 
scattering cross-section for the monopole on the 3K and 
diffuse photon backgrounds is just classical Thomson, 
valid oven for strong coupling: (Tt = 87ra^/3M^ ~ 
2 X iO-'*^ (M/10i°GeV)"2 ^.^^2. The resulting mean free 
path for inverse Compton scattering is many orders of 
magnitude larger than the Hubble size of the Universe. 

The relativistic monopole flux is simply 

Fm = cum/^tt ~ 2 X 10-^^ {M /10^°GeV f{lH/^cf (9) 

per cm^-s-sr, which compares favorably with the inte- 
grated flux above 10^" cV, Jdata(> lO^^eV) - 2 x 
10~^''cm~^s~^sr~^, and is comfortably below Parker's 
upper bound Fparker = 10~^^cm~^s~^sr~^ for a cosmic 
monopole flux. 

Signatures for EECR monopoles in our atmosphere 
and in ice are discussed at length in [29]. Of particular 
interest as a model for the super-GZK primaries is the 
"baryonic monopole" [30,29]. It is a bound state of three 
colored monopoles, confined by chromomagntic strings. 
Chromomagnetic strings can stretch when excited, but 
cannot break into quark- antiquark pairs. On first in- 
teraction in the atmosphere, the baryonic monopole of 
mass M and energy E = 7M stretches to create a huge 
geometrical cross-section of order 'j/Aq^y) ~ (l/lO®) ^ 
10^ mb. Consequently, nearly all of the initial monopole 
energy is transferred to the atmospheric shower in a very 
short distance. A recent simulation of the baryonic- 
monopole showed good agreement with the lateral muon 
and hadron content of the highest-energy Yakutsk event, 
but less than good agreement with the longitudinal pro- 
file of the Fly's Eye event [31]. 

Any confirmed directional pairing of events would ap- 
pear difficult to achieve with the monopole model. Also, 
in the context of a model of the Galactic magnetic field, 
it has been shown that some memory of the local spi- 
ral arm direction, and an energy spectrum flatter than 
the observed one, are expected in the data if monopoles 
are the primaries [32]. Although the directional crite- 
rion may not survive inclusion of cxtragalactic fields, the 
flatness of the spectrum probably does. 



C. Exotic physical laws 

The most remarkable proposals posit a breakdown of 
Lorentz invariance [33] or a breakdown of general rel- 
ativity above some high scale. In the case of broken 
GR, it could be spacetime fluctuations (expected in a 
theory of quantum gravity) wiggling the on-shell dis- 
persion relation [34] or appearing as 1/Mp operators 
which alter the physics. String theory provides moti- 
vation for GR-breaking at a possibly lower scale, the 
string scale M5. Certainly the energy window of EECRs 
(10^°eV /Mp ~ 10~^) is beyond that of terrestrial accel- 
erators, and so ripe for speculation on new high-energy 
physics. 

Signatures of models with no photo-pion production 

above Eqzk include the absence of a proton pile-up below 
Egzk, the absence of a cosmogenic neutrino flux, and 
possibly undeflected pointing of the primary back to its 
source. 



D. Neutrino primaries 

Turning to the possibility that the primaries may be 
neutrinos, one encounters an immediate obstacle: the SM 
neutrino cross-section is down from that of an electro- 
magnetic or hadronic interaction by six orders of magni- 
tude. This implies a low air-shower rate, and an accu- 
mulation of events at low altitudes ( "penetrating" events) 
where the target density is highest. On the other hand, 
the neutrino-primary hypothesis is supported by the ob- 
served clustering discussed earlier. Two solutions to the 
small cross-section problem for primary neutrinos have 
been proposed. 

1. Neutrino annihilation to Z-bursts 

Here it is proposed that the primary particles which 
propagate across cosmic distances above the GZK cutoff 
energy are neutrinos, which then annihilate with the cos- 
mic neutrino background (CNB) within the GZK zone 
{D < Dgzk ^ 50 Mpc) to create a "local" flux of nu- 
cleons and photons above Eqzk, as shown in Fig. (4). 
It was noted many years ago that a cosmic ray neutrino 
arriving at earth from a cosmically distant source has 
an annihilation probability on the relic-neutrino back- 
ground of roughly 3.0hQ^^% (neglecting cosmic expan- 
sion) [35]. The probability for a neutrino with resonant 
energy to annihilate to a .Z-burst within distance Dqzk 
is then 2.5 x 10~* for a homogeneous CNB [36]. The an- 
nihilation rate depends upon the CNB density, reliably 
predicted in the mean by Big Bang cosmology, and on 
the Standard Model (SM) of particle physics. The local 
annihilation rate is larger if our matter-rich portion of 
the Universe clusters neutrinos [36,37], or if there is an 
intrinsic CP- violating f — D asymmetry [38] . 
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Figure 4: Schematic diagram showing the production of a 
Z-burst resulting from the resonant annihilation of a cosmic- 
ray neutrino on a relic (anti)ncutrino. If the Z burst occurs 
within the GZK zone (~ 50 to 100 Mpc) and is directed to- 
wards the earth, then photons and nucleons with energy above 
the GZK cutoff may arrive at earth and initiate super-GZK 
air-showers. 

Each resonant neutrino annihilation produces a Z bo- 
son with a 70% branching ratio into hadrons known to 
include on average about one baryon-antibaryon pair, 
seventeen charged pions, and ten neutral pions [39] . The 
ten 7r°'s decay to produce twenty high-energy photons. 
For rr?,,y in the range ^ 0.1 to 2 cV. the energy in this reso- 
nant "Z-burst" is fortuitously situated sufficiently above 
£^GZK at 

= Mlllm^ = 4 (eV/m^) x lO^^eV (10) 

so as to produce photons and nucleons with energies ex- 
ceeding Eqzk* The mean energies of the 2 baryons 
and ^ 20 photons produced in the Z decay are easily esti- 
mated. Distributing the Z-burst energy among the mean 
multiplicity of 30 secondaries in Z-decay [39] , one has 

The photon energy is further reduced by an additional 
factor of 2 to account for their origin in two-body 7r° 
decay: 



*The resonant-energy width is narrow, reflecting the narrow 
width of the Z-boson: at FWHM AEr/Er ~ Fz/Mz = 3%. 



Even allowing for energy fluctuations about mean values, 
it is clear that in the Z-burst model the relevant neutrino 
mass cannot exceed ~ 2 eV. On the other hand, the neu- 
trino mass cannot be too light or the predicted primary 
energies will exceed the observed event energies. In this 
way, one obtains the approximate 0.1 eV lower limit on 
the neutrino mass, when allowance is made for an order 
of magnitude energy-loss for those secondaries traversing 
50 to 100 Mpc. 

The challenging issue of how experiments might actu- 
ally determine the absolute neutrino mass is discussed in 
[40]. 

If the Z-burst points in the direction of earth and oc- 
curs within the GZK distance, then one or more of the 
photons and nucleons in the burst may initiate a super- 
GZK air-shower at earth [36,37]. For a sufficient cos- 
mic neutrino flux, the hypothesis successfully explains 
the observed air showers above Eqzk- Comparisons of 
the model predictions to super-GZK data are available 
in [41]. 

The existence of neutrino mass in the desired range 
seems nearly guaranteed from the tritium decay upper 
bound [42] and the lower bounds inferred from the ter- 
restrial neutrino oscillation experiments. The simplest 
explanation for the atmospheric neutrino results is neu- 
trino oscillations driven by a mass squared difference of 
Jm^tm ^ 3 X 10~^eV^ [43], which implies a neutrino 
mass of at least 0.05 eV. Also, the recent LSND mea- 
surement appears to indicate a mass -squared difference 

(5mLSND ~ 0.2eV^ [44], from which one deduces a neu- 
trino mass of at least 0.5 eV. From these lower bounds 
on neutrino mass, one gets upper boimds on the Z burst 
energy of 10^'^ and 10^^ eV, respectively, just right for ex- 
tending the air-shower spectrum an order of magnitude 
or two beyond the GZK cutoff! 

A considerable cosmic neutrino flux above Eqzk is 
required for the Z-burst hypothesis to successfully ex- 
plain the super GZK events. The requirement is that 
the product of the resonant energy times the neutrino 
flux at the resonant energy per flavor, times the anni- 
hilation probability within the GZK zone (which may 
be as large as 0.025% to 1% due to neutrino cluster- 
ing) , times the; photon and nucleon multiplicity per burst 
(~ 20), is comparable to the observed flux at 10^° eV. 
The resulting requirement on the neutrino flux is roughly 
ErFi,.{Er) ^ 10"i**-^^Vcm^/s/sr. Such a neutrino flux 
at £" ~ 10^^ eV is directly measureable in a teraton 
(10^2 ton) detector like EUSO/OWL/AirWatch, and pos- 
sibly in a search for radio pulses produced by high energy 
neutrinos penetrating a small column density of matter 
in the limb of the moon [45] . 

While certainly large, this required neutrino flux vio- 
lates no existing limits. It has been pointed out [46] that 
this flux cannot extrapolate as E~^ to 10^^ eV, for then 
it would violate the Fly's Eye bound arising from nonob- 
servation of penetrating horizonal (i.e., neutrino induced) 
air showers at that energy. It has also been pointed out 
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[47] that local neutrino clustering is required to avoid 
generating a density of 30 MeV to 100 GeV photons 
in excess of the EGRET experimental bound [48], from 
distant Z-bursts undergoing electromagnetic cascading. 
The extreme-energy neutrino flux implied by the Z-burst 
model probably requires unusual source dynamics [49]. 
Among the reasons to hope that Nature obliges is that 
resonant neutrino annihilation provides the best hope at 
present to actually measure the relic neutrino density. 



2. Strong v Cross-section at ~ Eqzk 

It is interesting and suggestive that the observed 
EECR flux beyond Eqzk is well matched by the flux 
predicted for cosmogenic neutrinos. This is not a com- 
plete coincidence. With the GZK cutoff, any continued 
nucleon flux beyond Eqzk is degraded in energy, photo- 
producing pions which in turn decay to produce cosmo- 
genic neutrinos. The number of produced neutrinos com- 
pensates for their lesser energy, with the result that the 
neutrino flux matches well to the observed super-GZK 
flux. One may entertain the notion that the cosmogenic 
neutrinos are the super-GZK primaries, and that these 
neutrinos acquire a strong cross-section at ~ 10^° eV. 

Limits on the strength of the neutrino cross-section at 
10^" eV can be inferred from existing data. Heuristically, 
one argues as follows. The GZK process ensures that 
there is a flux of cosmogenic neutrinos at 10^'^ eV, with 
an easily calculated flux. If the neutrino cross-section 
were weak, an experiment looking for penetrating air- 
showers initiated by the cosmogenic neutrinos would see 
nothing. If the cross-section were strong enough, the neu- 
trino could not penetrate the atmosphere at all. So the 
fact that the Fly's Eye experiment saw no penetrating 
showers tells us that the neutrino cross-section is either 
strong or weak; the mid-range is excluded. The vertical 
depth of our atmosphere is ~ 1033 g/cm , and the 
horizontal depth Xh is about 36 times greater. In terms 
of the mean free path A of a particle with cross-section tr, 
one has x^jX = (7/1.6mb, and .Th/A = (T/44/ib. Thus, an 
estimate of the excluded cross-section is ^ 40/ib to Imb. 
A more careful calculation has been performed, with the 
result that ^ 20/ib to Imb is excluded [50]. Hypotheti- 
cal high-energy neutrino cross-sections in excess of a mb 
remain viable. 

The idea that neutrinos, indeed, all particles, may have 
a strong interaction at a high but observable energy scale 
is not new [51]. However, some recent ideas concerning 
new interactions relate well to the 10^" eV scale. One 
idea is that leptons are bound states of dual QCD gluons, 
which reveal themselves just above the electroweak (EW) 
scale at parton-parton -^s ~ TeV [52]. Another idea is 
that grand unification occurs precociously at \fs ^ TeV, 
because of extra dimensions or other reasons, and a neu- 
trino above this threshold becomes strongly-interacting 
via leptoquark resonances [53]. A third idea is that the 



exchange of a towers of Kaluza-Klein (KK) modes from 
extra compactifled dimensions lead to a strong neutrino 
cross-section above ^/s ^ TeV [54] . In all three cases, it 
is the combination of a low ^ TeV scale for radically new 
physics and a quickly rising spectrum of new states (pos- 
sibly increasing exponentially, p ^ eV'*/^o) that provides 
a rapid turn-on of a strong cross-section for the neutrino. 
Through unitarity, the new threshold at ^ TeV has con- 
sequences for cross-sections at lower energies [55], but 
they are not dramatic. 

The KK exchange model may fail [56] in that the KK 
modes couple to neutral currents, and the scattered neu- 
trino carries away 90% of the incident energy per in- 
teraction, thereby elongating the shower profile. But 
if the neutrino cross-section can be made large enough, 
^ 20 mb, then multiple scattering within a nucleus will 
effect a sufficiently large energy transfer and save the 
model [57]. Independent of the neutral current issue, 
the dual QCD and TeV-scale unification models seem 
to provide viable explanations of the super-GZK data. 
However, a recent calculation for the rate of rise of the 
low-scale unification cross-section in a string context is 
not encouraging [58]. 

Signatures for these models include directional point- 
ing back to the EECR source, longitudinal shower profiles 
differing somewhat from those of a proton or a gamma, 
and a strong correlation between observed energies and 
zenith angle. The latter signature should show an in- 
verse proportionality between the neutrino-air interac- 
tion length and the rising (with energy) neutrino cross- 
section. 



IV. MODEL SIGNATURES 

There are several telltale discriminators to be sought 
in higher statistics data. These will eventually eliminate 
most (perhaps all!) of the models so far proposed for the 
super-GZK events. We list some discriminating signa- 
tures and discuss them. 



A. Small-scale anisotropies and pointing 

The discriminatory power of small-scale clustering was 
discussed already in §11 A. Here we add some detail to 
the discussion. In traversing a distance D, a charged par- 
ticle interacting with magnetic domains having coherence 
length A will bend through an energy-dependent angle^ 



'On average, half of the interactions of a super-GZK nucleon 
with the CMB change the isospin. At energies for which cr 
of the neutron is small compared to the interaction mfp of 
~ 6 Mpc, the neutron decays back to a proton with negligi- 
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5e ~ 0.5° X :?^^7^Mpc Ampc . (13) 

-C'20 

Here BnG is the magnetic field in units of nanogauss, 
E20 and Z are the particle energy in units of 10^° eV 
and charge, and the lengths D and A are given in units 
of Mpc. It is thought likely that coherent extragalac- 
tic fields are nanogauss in magnitude [59], in which case 
super-GZK primaries from ~ 50 Mpc will typically bend 
only a few degrees (but note that protons at 10^^ eV will 
bend through ~ 30°). Thus, local models either postu- 
late many invisible sources isotropically-distributed with 
respect to the Galaxy to provide the roughly isotropic 
flux observed above Eqzk, or postulate a large extra- 
galactic magnetic field to isotropizc over our Northern 
Hemisphere the highest-energy particles from a small 
number of sources [60]. Among the latter category, some 
models postulate helium or iron nuclei as the primaries, 
to increase the bending by the charge factors 2 and 26, 
respectively. For those models invoking randomly dis- 
tributed, decaying super-massive particles (SMPs) [61] 
or topological defects (TDs) [62] as sources, and models 
invoking a large magnetic field with considerable inco- 
herent component, one may expect a nearly chance dis- 
tribution of observed events on the sky. However, there 
may be some clustering even in these models, due to pos- 
sible small-scale density fluctuations in the local SMP or 
TD distributions [63], or due to possible caustics in the 
projection of large-scale extragalactic magnetic fields on 
our sky (assuming the incoherent magnetic fields are suf- 
ficiently small) [64]. In the SMP and TD models, a high 
photon fraction in the primary composition further en- 
hances clustering possibilities. 

From the point of view of opening a new window to 
astronomy, those models in which the primaries do point 
back to their distant, active sources are the most inter- 
esting. These models are few in number. They arc the 
Z-burst model, the strongly-interacting neutrino model, 
and the quantum- gravity/LI- violating models. 

B. Large-scale anisotropics 

On large scales, one seeks associations of the CR direc- 
tions with the Galactic halo (to be revealed by a dipole 
anisotropy favoring the direction of the Galactic cen- 
ter) or the local galactic magnetic field, with matter 
distributions in nearby galactic or super-galactic clus- 
ters such Virgo, or with possible large coherent galac- 
tic or extragalactic magnetic fields [65]. For large-scale 



blc energy loss and the bending-angle formula is unchanged. 
However, at the energy 6 x 10^°, ct for the neutron is compa- 
rable to the interaction mfp, so at higher energies the nucleon 
bending-angle is reduced by ~ 2. 



studies, the Southern hemisphere Auger experiment will 
prove invaluable for several reasons. It offers coverage 
of potential sources and matter distributions, and galac- 
tic and extragalactic magnetic fields, not available from 
the North. Moreover, it offers a view of our Galactic 
center which will provide a North-South dipole discrimi- 
nator for or against a halo-centered population of sources 
such as magnetars or halo-bound SMPs. Southern Auger 
will also discriminate the M87-source model [23] wherein 
EECRs are channeled by a hypothesized galactic mag- 
netic wind into the Northern hemisphere, and the Ccn-A 
source model [24] which also yields a dipole anisotropy. 
Of course, an orbiting experiment with 47r vision like 
EUSO/OWL will be an even better instrument for mul- 
tipole analysis. 

C. Energy-direction-time correlations 

Because bending of charged particle trajectories by in- 
tervening magnetic fields increases as particle energy de- 
creases, one may learn about the strength and geometry 
of extragalactic fields from relative time delays and angu- 
lar correlations of particles from a common source. One 
may also learn about the source. Quantitatively, the in- 
crease in path length due to bending leads to a relative 
increase in travel time of 5t/t ~ (<5^)^, for small bend- 
ing angle. Adding the contributions from the coherent 
magnetic domains then yields 

5i~ 300DMpc f^l^) yrs (14) 

V -^20 / 

for the time delay. The time separation at earth is ob- 
tained by taking differences in eq. (14); to first order in 
SE it is already large: 

, , «nnn ( ZB^V ( 5E\ 

- i2 ~ 600Dmpc — ^ y^s. (15) 

V -'^20 / V-^20/ 

The correlation in energy and time becomes even more 
significant when it is remembered that the higher energy 
primary has an even higher mean energy in transit, before 
losses on the 2.7K background. 

Surprisingly, one AGASA event-pair has the higher- 
energy 1.06 X 10^" primary arriving about 3 years after a 
0.44 X 10^° primary. Assuming these primaries originate 
form a common source, a possible explanation is that 
their source has a duration of at least 3/(1 -I- z) years 
(the red-shift factor is due to time-dilation). Such a 
long-lived source does not occur in one-time burst models 
(e.g. GRBs). It may occur in decay models (e.g. SMPs. 
TDs, monopolonium) if sub-clustering exists on small- 
angular scales withing halo clusters [66]. A recent paper 
[67] notes that such sub-clusters may actually be observ- 
able as micro-lenses for stars and background galaxies. 
Counter-intiutive reverse pairing with the earlier arrival 
time for the lower energy charged-primary can also occur 
in certain magnetic field configurations, as shown in [68] . 
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D. Composition of the primaries 

Another signature to be sought is a statistical iden- 
tification of the nature of the primaries as a function of 
their energies. Three methods have been identified to dis- 
tinguish photon-initiated showers from hadron-initiated 
showers. One method rehes on the longitudinal profile 
of the event, particularly the depth at shower maximum 
a^max- The longitudinal profile of the Fly's Eye event at 
3 X 10^" cV (well-measured by its nitrogen fluorescence 
trail) is ill-fit by a photon primary, well-fit by an iron pri- 
mary, and somewhat fit by a proton primary [69]. The 
second method relies on measurement of muon number, 
with a high muon number purportedly favoring a hadron 
primary over a photon, and vice versa. A recent study 
of the muon content of showers above 10"'^^ eV seems to 
favor nucleons over photons [70]. However, caution is 
warranted with this method, in that some simulations 
show little difference in the muon-content of showers from 
nucleon primaries versus photon primaries [71] ; and pho- 
tons have themselves a significant partonic component at 
high energy. Nevertheless, it is safe to say that photon 
primaries appear disfavored at present, but more data is 
needed before conclusions are drawn. 

The third method of gamma identification will will 
rely on a predicted characteristic N-S vs. E-W gamma 
asymmetry. This quadrupole asymmetry results from the 
polarization-dependent interaction of the gammas with 
the earth's magnetic field. 

It is interesting to mention that the highest-energy 
Yakutsk event has an unusually high muon number; the 
only model so far which successfully explains this invokes 
a magnetic baryonic-monopole as the primary particle 
[31]. Basically, a relativistic baryonic-monopole of mass 
M and incident energy E showers like a giant nucleus of 
atomic number A^g ~ M/Aqcd ~ 10'^ (M/PeV) and the 
same 7 = E/M. Intermediate mass monopoles there- 
fore generate many, many charged pions which decay to 
muons. 

All models wherein the primaries arise from QCD jets 
produce many more pions than baryons. The neutral 
pions in turn produce gammas. In these models, the 
ratio of gammas to baryons is typically of order ten 
at the source. Even allowing for the shorter attenua- 
tion length of photons relative to nucleons, the measure- 
ment of the primary compositon for super-GZK events 
becomes an excellent discriminator for models with jet- 
producing sources. These models include Z-bursts, and 
decaying SMPs and TDs. 

E. Possible Smax energy cutofF 

The predicted of a cutoff at -Egzk is wrong. Does 
Nature provide an alternative cutoff within our reach? Or 
do the data continue beyond our reach? The shock-jock 



experts claim that it is difficult for conventional shock- 
acceleration mechanisms to produce ZeV proton energies; 
for this class of model, an i?max below 10^^ eV is certainly 
expected [72] . Decaying SMP models also have a natural 
cutoff, at half of the SMP mass. This could be as high 
as -Bmax ~ 10^^ eV for a long-lived GUT-mass particle, 
but would be lower for other postulated SMPs. In the 
Z-burst model there is a natural cutoff related to the tiny 
mass of neutrinos: .Bmax = A{eV/mv) ZeV. Implications 
from the atmospheric SK data are that this cutoff is at 
most 7 X 10^2 gv. 



F. CR flux above vs. below Eqzk 

A "smoothness" variable such as Rj = Fj{E > 
EG7.K)/FjiE < Eqzk) for each primary species 
j=nucleon, photon, iron nuclei, neutrino, etc. may be 
revealing [73]. For primaries with a GZK cutoff, F{E < 
Egzk) samples sources from the whole volume of the 
Universe, and may even include cascade products from 
F{E > -Egzk), whereas F{E > Egzk) samples just 
the GZK volume; for primaries without a GZK cutoff, 
F{E < Egzk) and F{E > Egzk) sample sources from 
the whole volume of the Universe. Lumps, bumps, and 
gaps in the spectrum near Egzk are a consequence of 
some models. For example, hadron and neutrino pile-ups 
just below Egzk are expected from the photo-pion pro- 
duction process which occurs above Egzk [74]. Other 
sources of neutrino pile-ups have also been suggested 
[36,75]. Present data show continuity. Smoothness stud- 
ies of various models require simulation, and are just be- 
ginning. 



G. Spectral index above Eqzk 

One means of achieving more events above Egzk is to 
postulate a flattening of the primary proton spectrum at 
highest energies. With more data, the extreme-energy 
spectrum will be measured. A flattening of the E^"^-"^ 
power law inferred from just below Egzk would indicate 
and constrain new sources. 



H. Measureable neutrino flux above Eqzk 

If there is a new source of primary neutrinos above 
Egzk, or if the neutrino cross-section becomes strong 
at super-GZK energies, then there is the possibility 
that the primary neutrino flux can be measured in an 
EUSO/OWL/AW-sized detector. This possiblity was 
discussed in section HID. 
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I. Diffuse ~GeV gamma-ray flux 

An upper limit on the diffuse gamma-ray flux between 
30 MeV and 100 GeV has been published by the EGRET 
experiment of the now defunct Gamma Ray Observa- 
tory [48] . This limit has serious discriminatory power for 
models where SMPs or TDs or extremely boosted mas- 
sive particles decay to quark-antiquark jets which then 
hadronization to produce the EECRs [76]. This is be- 
cause QCD jets via tt*^ production and decay produce 
very high energy gammas which initiate an electromag- 
netic cascade on the cosmic radio, microwave, infrared, 
and magnetic field backgrounds. For models with jet- 
production distributed over cosmic distances, such as 
some TD models and the Z-burst model with a homo- 
geneous distribution of relic neutrino targets, the cas- 
cade has the distance to reach completion; the end re- 
sult is gamma power (energy/time) in the EGRET range 
roughly an order of magnitude below the total power 
of the original sources, and comparable to the power in 
EECR neutrinos. Such models arc disfavored. More lo- 
cal models, such as SMPs bound to our halo, and the 
Z-burst model with a local over-density of relic neutri- 
nos, are not impacted at present, but may be tested with 
GLAST, the next generation gamma-ray observatory. 

V. CONCLUSIONS 

The ultimate explanation for the puzzles in EECRs 
will provide a surprise at a minimum, and possibly rad- 
ically new physics at a maximum. The extreme energies 
of events already observed cannot be approached by ter- 
restrial accelerators. Thus, there is ample motivation to 
build the next generation of CR detectors. Auger and 
Telescope Array, and to plan even farther beyond for a 
teraton detector like EUSO, OWL and Air Watch. At 
10^" cV, AGASA provides about an event per year, and 
HiRes about an event per month. Auger and TA will see 
two such events per week, while EUSO/OWL/AW has 
the potential to collect such an event every two hours. As 
the data sample grows, statistical studies will reveal sig- 
natures that discriminate among the many galactic and 
extragalactic sources so far proposed to resolve the EEGR 
puzzles. 

It is quite possible that neutrino primaries are respon- 
sible for the EECRs. If so, it appears that the weakly- 
interacting neutrino either grows a very strong cross- 
section at 10^" cV, or it annhilates on the relic neutrinos 
left-over from the hot phase of the Big Bang. Another 
possibility is that free magnetic charges exist and are the 
EECR primaries; for magnetic monopoles, the Universe 
is transparent and cosmic magnetic fields provide a nat- 
ural acceleration mechanism. 

Clearly, we live in exciting EECR times, and we pos- 
sess the technology to prove this is so. The resolution of 
our puzzles is forthcoming, as on-going and future exper- 



iments will provide us with the statistics to discriminate 
among the many interesting models. 
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